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Seventy crossings of the Earths bow shock by the THEMIS satellites have been used to study 
thermal electron heating in coUisionless, quasi-perpendicular shocks. It was found that the temper- 
ature increase of thermal electrons differed from the magnetic field increase by factors as great as 
three, that the parallel electron temperature increase was not produced by parallel electric fields, 
and that the parallel and perpendicular electron temperature increases were the same on the aver- 
age. It was also found that the perpendicular and parallel electron heating occurred simultaneously 
so that the isotropization time is the same as the heating time. These results cannot be explained 
by energy transfer from waves to electrons or by the motion of magnetized electrons through the 
shock. Electric field fiuctuations on the scale of the electron gyro-diameter were found to be of 
finite amplitude in the shock ramp, which requires that the electron trajectories be more random 
and chaotic than orderly and adiabatic. The data may be explained by the large amplitude electric 
field fiuctuations that demagnetize electrons as they move through the cross-shock electric field. 



INTRODUCTION 

During the half-century since the prediction [l|, @| 
and measurement [3| of coUisionless shocks, considerable 
study on mechanisms for converting supersonic plasma 
flow into thermal energy of plasma, fields, waves and ac- 
celerated particles has been expended. The removal of 
flow energy through dissipation, with the addition of dis- 
persion, is necessary to limit the nonlinear steepening of 
the shock HQ. 

The heating mechanisms for ions and electrons are be- 
lieved to be very different. While ion heating is fairly 
well understood, electron heating remains controversial. 
Proposed mechanisms for electron heating include (see 
e.g. id): 



• Adiabatic heating of magnetized electrons. The 
source of heating is the compression of the mag- 
netic field across the shock ramp. This reversible 
mechanism predicts heating perpendicular to the 
magnetic field {Te±). 

• Kinematic heating of magnetized electrons, where 
the cross-shock potential in the ramp gives a cor- 
rection to the adiabatic moment. This scenario is 
also dependent on magnetic compression, and leads 
to super- adiabatic perpendicular heating (Tej_). 

• Parallel heating (Te||) of magnetized electrons by 
the parallel electric field in the shock ramp. 

• Energy and momentum transfer from waves to elec- 
trons through processes such as anomalous resistiv- 
ity. Here the source is unstable currents or parti- 
cle distributions in the shock layer giving rise to 
e.g. lower-hybrid, ion-acoustic or whistler micro- 
instabilities and turbulence. Wave-particle inter- 
action causes momentum exchange, where waves 
gives energy to the particles (anomalous resistive 
heating). 



Of these mechanisms the first three are due to macro- 
scopic fields in the shock, while the latter is microscopic 
in nature. 

Through examination of 70 terrestrial bow shock cross- 
ings made by the THEMIS satellites, it is shown in this 
paper that none of these mechanisms adequately explain 
the electron heating observations. Instead, it will be 
shown that electron heating is achieved by turbulent 
transport of unmagnetized electrons through the cross- 
shock electric field. 

The paper is organized as follows. First the instru- 
mentation and data sets are presented. Then, the above- 
listed heating mechanisms are investigated, one by one, 
and it is shown that none of them are consistent with 
the experimental data. Finally, the model of electric 
field turbulence that demagnetizes electrons as they move 
through the cross-shock electric field is found to be ade- 
quate to explain the data. 



INSTRUMENTATION AND DATA SETS 

We use data from the fiuxgate magnetometer [l^ , the 
electron electrostatic analyzer 17], and the electric field 
instrument [l8| on the five THEMIS spacecraft to study 
bow shock crossings at times when the data transmission 
rate was 128 or 8192 samples per second and, usually, 
when three of the spacecraft encountered the Earth's bow 
shock within a short time interval. We selected quasi- 
perpendicular shocks, all of which were super-critical 
with respect to the Alfvenic Mach number. Ma- 

Electron temperatures are conventionally obtained 
from moments of the electron distribution function. 
In this study we, instead, use the electron tempera- 
ture estimated from fitting Maxwell-Boltzmann curves 
to the electron distributions, as illustrated in Figure [H 
In this Figure, the blue diamonds are perpendicu- 
lar (75° to 105°) pitch angle electrons measured dur- 



2 



2; 
o 

H 

O 



o 

H 



10" 



13 



10' 



m Pi 

I— ( 

Pi 

H 



10^ 



10" 



"I 1 

THEMIS C, June 26, 2008 1640 UT 



1 3 e V Maxwellian 




♦ ♦ 

l_ 



10 



100 

ENERGY, eV 



1000 



FIG. 1: Distribution function of 90° pitch angle electrons 
measured during a one second interval in the solar wind, 
and Maxwell-Boltzmann curves for temperatures of 7 eV and 
13 eV (13 eV is the moments estimate of the temperature for 
this data). 



ing a one-second-interval upstream of the Earth's bow 
shock. This distribution consists of a low energy ther- 
mal core and a supra-thermal power law of halo elec- 
trons. The Maxwell-Boltzmann fit to the thermal core is 
the red curve characterized by a 7 eV temperature. The 
moments-temperature of 13 eV is illustrated by the green 
Maxwellian. The moments-temperature is greater than 
the Maxwellian-temperature because the moments fit in- 
cludes contributions from the supra-thermal electrons. In 
the current data set, the moments-temperature exceeds 
the Maxwellian temperature by factors between about 
1.2 and 2.0. An additional reason for using Maxwellian- 
temperatures in this study is that they may be com- 
puted at a cadence much less than one second, while the 
moments-temperature is available only once every spin 
period (four seconds for the THEMIS spacecraft). 



HEATING BY ENERGY TRANSFER FROM 
WAVES TO THE ELECTRONS DOES NOT 
EXPLAIN THE DATA 



The idea that anomalous resistivity can be a channel 
for heating of plasma in shocks has been studied since 
the early days of collisionless shock research. In this sce- 
nario unstable currents in the shock transition layer, as 
well as unstable particle distributions and gradients in 
density, temperature and magnetic field, can drive insta- 
bilities which excite a spectrum of waves. These waves 
would then, through wave-particle interaction, lead to 
momentum and energy exchange between waves and par- 



ticles. While this process is microscopic (kinetic) the 
macroscopic impact is equivalent to a form of resistivity 
with an associated collision frequency. Since there are 
no binary collisions taking place, the resulting resistiv- 
ity is called " anomalous" . A formal derivation using the 
Vlasov equation led to an expression for the effective col- 
lision frequency 19 L A similar derivation from a gener- 
alized Ohm's law 20, 21 [, leads to the expression for the 
part of the cross-shock electric field due to electrostatic 
anomalous drag or resistivity 



= — 



< SnSEx > 
<n> 



(1) 



where n is the plasma density and is the electric field. 
Dj. provides the dissipation required for converting elec- 
tromagnetic energy to plasma energy and it has been 
used to show that wave heating in magnetic field recon- 
nection is unimportant 21]. 

In Equation [Jl the numerator includes the fluctuations 
of the electric field and density and their correlation. It is 
computed for a bow shock crossing on October 24, 2011 
in Figure [5] (this crossing is discussed in greater detail in 
figures and [7]) . Panel[2^ of this figure gives the 0.25 s 
average of the total magnetic field during the crossing, 
panel [2j) gives the 0.25 s average of the cross-shock 
electric field, and panel[2l; gives the 0.25 s average of Dx- 
Because Dx ^ Ex, wave heating of electrons across the 
bow shock is unimportant. Although the fiuctucations 
in the electric field are large compared to the mean field 
(see Figure [Tji), and the density fluctuations are large 
compared to the mean density (see Figure [5)d), Dx is 
small because the electric field and density fluctuations 
are uncor related. 



PERPENDICULAR HEATING OF MAGNETIZED 
ELECTRONS DOES NOT EXPLAIN THE DATA 



Early work by Goodrich and Scudder [6|, |]J-UJ, 123 
considered magnetized electrons. To a first approxima- 
tion, magnetized electrons conserve the fluid first adi- 
abatic invariant, Te±/B, across the shock ramp where 
Te± is the perpendicular Maxwellian-temperature and B 
is the magnetic field magnitude. 

The adiabaticity is defined as 

adiabaticity = {Te±^sh/Bsh)/{Te±,sw/ Bsw) (2) 

where sw and sh refer to the upstream (solar wind) and 
downstream (sheath) parameters. For magnetized elec- 
trons the adiabaticity should be equal to one. For 70 
quasi-perpendicular {0Bn between 46° and 89.5°) shocks, 
the adiabaticity is plotted in Figure [3^ versus the Alfven 
Mach number. Ma- At individual crossings, thermal elec- 
trons may be heated more or less than adiabatically by 
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FIG. 2: Measurements across a bow shock on October 24, 
2011. Panel 2a gives the magnetic field magnitude, panel 2b 
presents the cross-shock electric field, and panel 2c gives the 
anomalous drag, D. Because D is much less than the cross- 
shock electric field, conversion of wave energy to particle heat- 
ing is relatively unimportant. 
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FIG. 3: Scatter plots of the adiabaticity (panel 3a) and ratio 
of parallel to perpendicular electron heating (panel 3b) for 70 
bow shock crossings. Panel 3a shows that the perpendicular 
electron heating was as much as three times larger or smaller 
than that expected for adiabatic compression. Panel 3b shows 
that the average parallel and perpendicular electron heating 
were the same. 



factors as great as three (although the average adiabatic- 
ity is 1.06 ± 0.52). Thus, the model of conservation of 
the first adiabatic invariant used to explain heating in the 
compressed magnetic field is inconsistent with this exper- 
imental data. The adiabaticity obtained using moments- 
temperatures has the same characteristic as does the data 
of Figure dSJa. 



ELECTRON HEATING IN PARALLEL 
ELECTRIC FIELDS DOES NOT EXPLAIN THE 
DATA 

Parallel heating of magnetized electrons cannot be 
achieved by adiabatic compression. Instead we con- 
sider parallel heating by the parallel electric field. In 
the deHoffman- Teller (HT) frame [l^l, the perpendicu- 
lar electric field in the solar wind is zero, so the plasma 
moves parallel to the magnetic field. In an idealization, 
the perpendicular electric field is also zero through the 
shock, and downstream. In this case, the parallel heating 
can only be achieved by the parallel electric field in the 
ramp. To test this idealized model, a typical example of 
the electric field at a shock crossing is given in Figure 1^ . 
Panel 13^ gives the total magnetic field and the remain- 
ing panels give the HT electric field in a magnetic-field- 
aligned coordinate system in which panel |4ji is the com- 
ponent parallel to the instantaneous magnetic field. The 
electric field components in the HT-frame are all small 
upstream and downstream. The parallel component of 
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FIG. 4: The total magnetic field (panel 4a) across a bow 
shock and the three components of the electric field in the 
deHoffman- Teller frame, which is the frame in which the per- 
pendicular electric field is zero in the solar wind (at the right 
of the plot) . The fields are in a magnetic-field-aligned coor- 
dinate system with the electric field component parallel to B 
given in panel 4d. In this example, and for all of the 70 events, 
the parallel electric field in the shock ramp was small com- 
pared to the perpendicular field in this or any fixed coordinate 
system. 
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panel |4ji is zero to within the experimental uncertainties 
and, contrary to the idealized model, the perpendicular 
components in panels |3)d and are large in the ramp. 
It is noted that, for the event of Figure 21 the magnetic 
field was nearly in the spacecraft spin plane so the parallel 
electric field was largely measured by the long spin-plane 
electric field sensors and the large uncertainties associ- 
ated with the short spin-axis electric field sensor had only 
a small effect in the parallel field measurement. 

This data is representative of our observations that, ex- 
cept for occasional short duration bipolar parallel fields 
due to Debye scale electron holes, the parallel electric 
field in the shock ramp is much smaller than the per- 
pendicular field. Thus, the parallel electron heating does 
not arise from the parallel electric field in the ramp, but 
rather, from the effects of the perpendicular electric field. 

This fact is confirmed in Figure [3)d where it is shown 
that the parallel and perpendicular electron heating in 
each of the 70 events differed by at most 40% and was, 
on the average, equal. 

UNMAGNETIZED ELECTRONS 

Electrons would be unmagnetized if fluctuations of 
the electric fleld over the electron gyro-orbit were sig- 
nificant. To test if this is the case, fluctuations, S<j>, 
in the cross-shock potential (where (j> — J E^Vghdi and 
Vsh = 10 km/s) as well as density fluctuations, 6n, are 
used to compute 

Sn/n=[n(t + T) -n{t)]/ <n> (3) 

S^/T, = + r) - ^{t)]/ < Te > (4) 

where r is the time required for the shock to move a 
distance across the spacecraft equal to the electron gyro- 
diameter and the averages in the denominators are com- 
puted over the time interval from t to {t + t). The 
electrons are unmagnetized if Scji/T^ is of the order of 
one. Figure [5] presents data for a shock crossing by 
THEMIS E on October 24, 2011 (the same event as in 
Figure [2]). The field data was obtained with a sam- 
pling rate of 128 S/s, and the perpendicular temper- 
ature was measured every 0.25 s. The parameters for 
this shock were Ma = 5.5, dsn — 85°, the normal was 
"-dsi = (0.9991, —0.0235,0.0356) in despun spacecraft co- 
ordinates (essentially Geocentric Solar Ecliptic coordi- 
nates), and the local time of the crossing was 14:30 UT. 
Panel [5^ gives the total magnetic field, panel gives 
Sn/n and panel [5]: gives 6(j)/Te±. The quantity r in these 
computations is 70 msec, which is the time required for 
the shock to move the gyrodiameter of a 40 eV electron 
in a magnetic field of 50 nT. Through the shock ramp, 
the gyrodiameter changed by a factor of about two, so 
the computations of Sn/n and S(j)/Te were also done for 
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FIG. 5: The bow shock crossing that is also illustrated in Fig- 
ures [21 [6l and [71 This data is in the normal incidence frame 
(NIF) in which fluctuations over the electron gyro-diameter 
of the plasma density (panel ^jp) and the cross-shock elec- 
tric potential (panel Ht) are important in the shock ramp. 
The perpendicular electron heating (panel[5li) correlates with 
the electric field turbulence of panelist and the one-second- 
averaged electric field of panellSf. The parallel heating (panel 
15^), has the same amplitude and time dependence as does 
the perpendicular heating (within the 1-2 second time reso- 
lution of the parallel measurement). This indicates that the 
parallel and perpendicular heating result from the same phys- 
ical process. Their correlation with the potential fluctuations 
of panel ISj; and the electric field of panel 15}' shows that this 
process involves unmagnetized electrons moving in the cross- 
shock electric field. 

T values ranging over a factor of two with similar results 
to those presented in Figure l5] These results are that 
S(j)/Tej_ is of the order of one in several locations in the 
shock ramp, so the electron trajectories must be more 
random and chaotic than orderly and adiabatic, and the 
thermal heating of electrons must be of a stochastic na- 
ture. 

Data were collected at 128 samples/second through the 
entire crossing of FigurelH corresponding to observations 
of fluctuations below 50 Hz. For a period of 0.3 s in the 
ramp, data were also collected at 8192 samples/second. 
The power spectrum of the cross-shock electric field for 
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FIG. 6: Spectrum of the electric field fluctuations in the bow 
shock crossing of Figure (5] Although there is power at and 
below the lower hybrid frequency, the dominant wave power 
is in the whistler frequency range. 



this time interval is given in Figure |6l Fluctuations near 
the lower hybrid frequency and in the whistler mode 
range were observed, with the average waveform ampli- 
tude in the whistler mode range being a factor of 2.6 
greater than the average amplitude covered by the 128 
Hz sampling rate. Thus, the 128 Hz fluctuations, were 
multiplied by 3.6 to produce the plot of Figure[5t. 

These fluctuations in 5(j)/Tf,±_ are not bipolar electrostatic 
structures [24| at Debye scales in the ramp. Such stuc- 
tures are observed occasionally in this and other shock 
crossings, but their integrated potentials are small 25 1. 

The data of Figure [5] are in the Normal Incidence 
Frame (NIF). The NIF is defined as a shock stationary 
frame where the incident plasma flow is aligned along the 
shock normal. The fields in this frame are given in Fig- 
ure[7l in which the top three panels are the components of 
B and the bottom three panels are the components of E. 
The x-direction is normal to the shock. The co-planarity 
plane includes Bz and Ex- In the solar wind, to the left 
of the plot. Ex ^ (panel [TH), while in the ramp and 
later it is clearly positive. Thus, there is a cross-shock 
electric field and a net potential. The combination of fi- 
nite amplitude electric field fluctuations and associated 
potential in the presence of the DC cross-shock electric 
field provides the mechanism for heating the electrons as 
the following discussion illustrates. 

Figure [5ji gives the perpendicular Maxwellian- 
temperature, Figure [5]3 gives the parallel Maxwellian- 
temperature, and Figure [Sf gives the 1 s averaged cross- 
shock electric field. The correlation between the perpen- 
dicular heating (panel [SJi), the electric field fluctuations 
(panel [SJ:) and the cross-shock DC electric field (panel 
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FIG. 7: The three components of the electric and magnetic 
field during the bow shock crossing illustrated in the previous 
figures. The data are in the normal incidence frame (NIF) 
such that the solar wind plasma flow is in the direction of the 
shock normal, X. The magnetic field component of panel [Tj; 
and the normal, determined from a minimum variance analy- 
sis, are in the co-planarity plane, and the cross-shock electric 
field (panel d) is small upstream but positive in the ramp and 
downstream. 



[Sf) is remarkable. This shows that strong electric field 
fluctuations of finite order, in association with the non- 
zero cross-shock electric field causes the electron heating 
in the shock ramp. In the latter portion of the plot, 
there continues to be a positive electric field (panel [5f ) 
although the electron temperatures (panels [SJi and [5^) 
do not increase. Thus, in the absence of turbulence, the 
perpendicular electric field caused nothing more than an 
E X B/B^ drift. The parallel electric field, being small 
or zero, also did not heat electrons at this time. 

The parallel electron temperature, obtained for elec- 
trons with pitch angles between 0° and 30°, is shown 
in Figure [5^. Within the 1-2 second time resolution of 
this measurement, the parallel and perpendicular tem- 
peratures were equal and their rise times were the same. 
This indicates that the parallel and perpendicular heat- 
ing result from the same physical process. Their corre- 
lation with the potential fluctuations of panel [5j; shows 
that this process involves unmagnetized electrons in the 
cross-shock electric field of Figure [7ji. 

The electron perpendicular and parallel temperatures 
also correlate with the magnetic field of panel [5^. While 
this might suggest that the compression of the magnetic 
field plays a role in the electron heating, it is known 
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FIG. 8: Example of a bow shock crossing during which per- 
pendicular electrons were heated in the shock ramp (panel 
d) coincident with large electric potential fluctuations (panel 
c) and the magnetic field increase in the ramp (panel a). 
Slightly later, in the downstream direction, the magnetic field 
decreased by a factor of five (panel a) in the absence of large 
amplitude potential fluctuations (panel c). At this time the 
perpendicular electron temperature (panel d) remained un- 
changed. 



from the adiabaticity plot of Figure [3] that this is not the 
case. Another example of the relative unimportance of 
the magnetic field for electron heating is given in Figure [S] 
in which the temperature increase (panel |5Jl) correlates 
with the electric field turbulence (panel [S}:) and also with 
the magnetic field increase (panel |8^). However, slightly 
later, the magnetic field decreases by a factor of five and 
there is no change in the perpendicular electron temper- 
ature (panel IHJi), so the magnetic field magnitude is not 
controlling the electron perpendicular temperature. The 
adiabaticity for this event is 2.24. It is also noted that 
large values of dnjn and 5(j)/Te±_ have been reported for 
a different event 1261. 



DISCUSSION 

It is shown that neither anomalous resistivity nor mag- 
netized electron fiow are capable of explaining thermal 
electron heating at terrestrial bow shock crossings ob- 
served by the THEMIS satellite. Instead, the combina- 
tion of strong electric field turbulence and the cross-shock 
electric field in the NIF are required to understand why; 

• The perpendicular electron temperature increase 
differs from the magnetic field increase by factors 
as great as three. 

• The parallel temperature increase is much greater 



than can be provided by the parallel electric field. 

• The perpendicular and parallel temperature in- 
creases across the shock are, on the average, equal 
and they increases simultaneously. 

Of the 70 bow shock crossings in the study, 67 had 
some level of electric field fluctucations similar to those 
described above. At the rate of 128 samples/second, 
75% of the crossings had turbulence levels greater than 
10 mV/m. There was a weak dependence of the > 

10 mV/m turbulence on the solar wind electron tempera- 
ture and magnetic field, with large amplitude turbulence 
found in all six events having a solar wind temperature 
above 12 eV and 10 of 11 events having a solar wind 
magnetic field greater than 10 nT. 

Because the events were chosen selectively and be- 
cause the apogees of the THEMIS satellites were about 

11 Earth radii compared to a typical bow shock loca- 
tion of ~ 15 Earth radii, the analyzed events may not 
be typical of the bow shock physics for a randomly se- 
lected set of crossings over a wider range of altitudes. To 
test this possibility, bow shock crossings by the Cluster 
2 satellite during January, 2003, were examined for evi- 
dence of electric field turbulence. (Cluster 2 transmitted 
0-180 Hz bandwidth data at a lower data rate, so this 
aliased data was appropriate for observing whether tur- 
bulence was present, but not for quantitatively analyzing 
it.) The 140 clean Cluster bow shock crossings during 
this month occurred over radial distances of 12.9 to 18.9 
Earth radii, and 90% of them had fiuctuation amplitudes 
> 10 mV/m. Thus, the THEMIS events in the current 
study were not anomalous and the results of this study 
may be extrapolated to all terrestrial perpendicular bow 
shocks. 

We finally arrive at the following new scenario for ther- 
mal heating of electrons in super-critical shocks, consis- 
tent with the observations presented above. Instabilities 
due to gradients, currents and possibly unstable particle 
distributions in the shock ramp are driving strong non- 
linear fiuctuations, primarily in the electric field. When a 
magnetized electron drifts into regions of strong potential 
fiuctucations on the scale of a gyro orbit, 5(j)/Tf.± ^ 1, 
the electrons demagnetize in this local region, 6(j)/Te±, 
Te±, Tgii and the DC electric field correlate and most of 
the electron heating occurs. This scenario is required be- 
cause it has been shown that neither magnetized electron 
heating nor heating due to anomalous resistivity is con- 
sistent with the experimental data. It should however be 
pointed out that such strong potential fields are likely to 
trap electrons, which could lead to secondary effects not 
explored here. 

Because astrophysical shocks should be stronger than 
terrestrial bow shocks, the mechanism of electron heating 
by the combination of strong electric field turbulence and 
a cross-shock electric field may be applicable to all of 
astrophysics. 
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